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INTRODUCTION 

Increasing concerns about sustainable energy generation have triggered considerable advancement 
in the tidal stream industry. The tidal stream industry is in the intial stage of converting the full-scale tested 
prototypes to a commercial level, therefore the planning of tidal stream sites needs to follow a strategy 10 
‘maximize their potential [1]. Among the marine renewable energies, tidal stream, as well as tidal range 
fechnology, are the most promising renewable sources due to their highly predictable potential [2], [3]. 
Because of the increased prediction and accuracy rate of the tides, the output power of a particular tidal 
power plant is also highly predictable, which is an essential parameter in the cost-benefit evaluation phase 
Compared to wind energy generation, the tidal stream flow highly depends upon the phases of the moon, and 
sunt, and therefore tidal power prediction is comparatively more accurate [4]. Furthermore, the density of 
seawater is almost 800 times higher than the density of air, in this way the power generation by a Tidal 
energy converter (TEC) will be higher than the wind machine at an appropriate rated speed. Moreover, the 
extreme velocity of the wind can damage the wind machines, however, the TECs do not incorporate such 
damage because of the absence of extreme tidal flow velocities. Apart from the advantages of tidal stream 
technologies. they can lay minor negative impacts on the marine environment. The major drawbacks include 
the impact on the biota atthe plant site, morphodynamics, and sediment transport. It is therefore important 10 
hhave a detailed impact assessment before selecting the tidal steam energy site 
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2, LITERATURE REVIEW 

Oceans around the world possess huge potential to harness elean energy and take a great part in the 
decarbonization of energy systems. ‘The tidal power potential was theoretically computed to be approximately 
7800 TWhiyear worldwide [5]. A large segment of such a potential is yet to be explored. Iti estimated by 
the year 2080, a cumulative replacement CO2 emissions can be reduced by 1.0 billion tonnes thereby 
generating 337 GW of clean energy from waves and tides [3]. Marine energy utilization contributes to the 
sustainable energy supply of the world. In the early stage installed marine energy capacity, considering wave 
and tidal energy from the year 2000 to 2010, Europe was dominant with an installed capacity of 241 MW, 
covering 91% of the World's installed marine energy capacity, as shown in Figure | [6]. Whereas, share from 
Nom America and Asia was 20 MW (7.5%) and 4.5 MW (1:8), respectively. Afterward, there is no such 
remarkable progress in the installation of marine energy systems in Europe from the year 2010 t0 2016 
‘There was a cumulative addition of 6 MW only during that period, and hence decreasing its Worldwide share 
from 91% to 46% {6]. However, considerable growth in installed capacity can be seen in Asia, which 
Increased from 4.5 MW to 262 MW, being dominant to have 49% share in the world’s installed marine 
energy capacity [6]. Such a huge inerease in installed capacity share in Asia is mainly because the Sibwal 
Lake Tidal power plant was commissioned in the vicinity of Seol, South Korea in 2011. Afterward, the 
number of marine projects are being promoted and supported by the Korean government due to its 
geographical location to receive the abundance of tidal flow from its three sites. I is estimated that tidal 
current energy resources along the southwest coast can generate up to 1000 MW of clean power production 
7]. There has been continuous progress in research and innovation to bring advancements in Marine 
renewable energy, in parallel with other renewable energy sources. 

‘Technological advancement of marine energy systems is the key feature to bring it to the 
commercial level. Marine energy researchers have tied their best to bring innovation in marine energy 
technology at an appreciatable level when considering other renewable energy systems, which can be seen 
from the number of patents filed worldwide in Figure 2. The majority of the patents for renewable energy 
systems are solar-hased energy systems [6]. However, the maximum number of patents filed in marine 
energy was 1,533 in the year 2010 with an increase in relative change of 552%. The majority of the patents 
filed for renewable energy systems were from China, covering 66.6% of the total fied patents in the year 
2016 [6]. The majority of the countries abundant with tidal current and the ocean current potential include 
South Korea, China, Norway, France, Canada, Australia, New Zealand, India, the Russi, the UK, and the 
USA [8]. However, other countries such as Malaysia, Spain, Taiwan, Indonesia, and Iran also have tidal 
current energy sites from where a significant amount of tidal energy can be generated [5], [9|-11]. The 
considerable attention has been given for research in tidal energy exploration in these countries, and the new 
sites are explored for TEC installation and sea tral. 

‘China is dominant to have an abundance of tidal stream energy generation sources in Asia. Having 
some of the tidal energy sites to be observed with tidal current speed more than 3.06 mis while the majority 
inthe range 2.06 to 3.06 mis, China is estimated to generate about 61.3 TWh/year of energy [12]. 
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Figure 1. Marine energy capacity installed worldwide from the year 2000 to 2016 (61 
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Figure 2. The number of patents filed in the World for renewable energy systems from the year 2000 to 2016 [6], 


2A. 


‘dal energy prospect in the Northern Indian Ocean 

The Nothern Indian Ocean is divided into two seas namely the Arabian sea and the Bay of Bengal. It is 
shared by some of the South Asian countries including Myanmar, Bangladesh, India, Sri Lanka, Maldives. 
Pakistan, Iran, Oman, Yemen, and Somalia [13]. However, this section is focused to discuss the tidal energy 
prospect of Bangladesh, India, Pakistan, and Srilanka, as indicated by the boundary line in Figure 3. The Indian 
‘ocean is mainly distinguished due to the reversed monsoon season and ocean currents [14], [15]. The surface 
‘ocean currents prediction is shown in Figure 3, covering the portion of the Bangladesh and the Arabian Sea that 
is inthe vicinity of India. 


Figure 3. The surface current speed in knots in the Northern Indian Ocean is predicted on 14th June 2021 


From the appraisal in Figure 3, Sri Lanka is noted to be the most dominant to have the potential tidal 
‘current. The sea surface current speed ranges 4 to 5 knots during the different times of a neap tide along all of its 
coastline. Whereas the main sites for India include the sites close to Tamil Nadu, Machlipatnam, Brahmagiri, 
and Surat, in these locations, the surface current speed up to 5 knots was noted during the neap tide on 14th June 
2021. However, the potential sites can be further investigated by developing hydrodynamic models for these 
locations. In Pakistan, the tidal site near Karachi was observed to have a surface current speed of 3 knots. 

Prediction for the wave height is shown in Figure 4. The wave energy is dissipated slowly over a very 
Jong distance from its origin. The distribution of wave heights in Figure 4 also reflects the patterns of wind in 
the Norther Indian ocean, The energy dissipated by the waves undergoes daily as well as seasonal variations 
as itis originated from the wind [16]. It is clear from Figure 4, that the majority of the Indian coastline is 
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dominant with a wave height of around 3 m. Although the waves are originated due to wind, wave power is still 
‘more predictable than wind energy generation [17]. Wave energy density is measured in kW/m which stands for 
the energy density per wave erest. A typical wave energy converter is optimized for the energy density ranges 
from 15 to 35 kW/m, whereas 20 kw/m is considered as a feasible site [18]. However, in tidal current-based 
energy generation, the tidal current speed of Ito 1.5 ms can be regarded as a good potential site [19], 

‘A detailed hydrodynamic study can be carried out along the coastline of India to identify the exact, 
locations to be the potential tidal current energy sites. Tidal energy development also needs the support of 
govemment officials and policymakers. Most countries with remarkable progress in tidal energy generation 
have set targets for renewable and ocean energy generation, to promote tidal energy development 


Figure 4. The wave height in the Northern Indian Ocean is predicted on 14 June 2021 


3. HYDRODYNAMIC MODELLING FOR TIDAL ENERGY EXTRACTION 

Hydrodynamic modeling for a particular location is accomplished using numerical models suchlike 
POM, Mike 21, and TELEMAC. However, for tidal energy assessment, Del 3D is widely used by the 
research community to solve the hydrodynamic problem with high resolution (1), [9 {10}, [20}-(23]. The 
software is originally developed by Delft Hydraulics, and it is an open-source code freely available for 
everyone to use for academic studies. It solves the shallow water equations, by discretizing the model area 
that is covered by the curvilinear grid. For the hatter result, the constructed grid should be well structured and 
well orthogonalized. The flow is mainly described by the two main variables namely water level and 
velocity. The density/pressure points (i.e water level points) are assigned at the center of each grid cell, 
Whereas velocity is assigned perpendicular to grid sides. Such arrangement of the orthogonal physical 
quantities in a computational grid is called Arakawa C-grid (Delft3D). The time step condition to integrate 
the shallow water equation on the orthogonal grid ean be computed based on the Courant number, which is 
expressed in (1) [24] 


Flyase = 200VGH [Et <1 a 


In the above equation, the time step is represented by At, and H represents the total water depth, g denotes 
gravitational acceleration, and Ax and Ay stand for smallest grid distances [24] 

To solve the shallow water equations for a hydrodynamic problem, the incompressible fluid is 
assumed during the simulation. The continuity and momentum equations solved by Delfi3D-FLOW for a 
hydrodynamic problem are given below (Delft3D). 


Tndonesian J Elec Eng & Comp Sci, Vol. 29, No. 2, February 2023: 1201-1212 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 1205 


roa) 
de + 2 + yy TU, 
2p 4 4 yO 
det AB + TY 8) 
Mee 5 sate 4 EME pore Ade De+R ® 


According to (2) and (3) defines conservation of mass and momentum, respectively. According to 
(4) is used for salinity and transport. Where "f denotes water level, “d” is relative water depth, "U” is the 
velocity in X direction and “V" is the velocity in Y direction, “Q” is the mass intensity in single grid cell, “f" 
denotes Coriolis constant, “wy,” denotes horizontal eddy viscosity, “p” stands for seawater density, “t, 
shows surface wind stress, “t,” is the bottom shear stress, and "Ris force on single grid cell 

‘The tidal energy assessment is mainly accomplished by two methods, ie based on the actuator dise 
theory method or based on the momentum sink method [25]. Tidal energy assessment by Delft3D implies the 
‘momentum sink approach [22]. This is based on energy loss in the Navier Stokes equation when any obstacle 
is present in the water flow. In this way, an additional term for energy loss is added to the momentum 
equation. That obstacle (in Delft3D it is done by adding porous plate as an additional parameter) works in the 
same as an array of turbines will reduce the pressure of flow across the turbine (ie momentum loss or 
‘momentum sink). In this way, energy extraction from the flow can be modeled. By adding the energy loss 
coefficient to the momentum equation, momentum in (3) in x-direction will become as: 


ay yy ya 2 (1,20) 4 2 (20) _ 1008 

eeu FV = (03) +50) 3) +h Ma 2) 
as shown in (5), the additional term fx stands forthe horizontal Reynolds stress, and Myx shows momentum 
sink in x-direction that is perpendicular tothe incoming flow [22 Similarly, the momentum sink term can be 
audded in the y and z direction when performing 2D and 3D tidal energy extraction modeling 


4. CASE STUDY OF SARAWAK MALAYSIA 

‘Malaysian tidal power development still remains in early stage of feasibility investigation. Several 
research studies have been carried out in the lst decade to identify suitable sites as well as to investigate the 
energy extraction potential of some of the identified sites. The first research investigation to assess the 
capability of Malaysian tidal power was performed by (26). in which the three main sites, namely Pulau 
Jambongan, Sibu, and Kota Belud, were identified as the main locations having the great potential for tidal 
stream energy hamessing. Afterward, other researchers also contributed to assessing the potential of tidal 
sites by developing the hydrodynamic models using Delft3D numerical tools. [20] explored the tidal current 
energy sites and seabed roughness at Stat of Malacca, they proposed the Pulau Pangkor as the mos suitable 
location for installing TEC. Bonar er a (27] developed the depth-averaged hydrodynamic model for further 
assessing the tidal power generation capability a the Strait of Malacca, by considering a row of tidal turbines 
as the actuator dse. Later on, [21] developed a Delft3D hydrodynamic model forthe Malacca Strait focusing 
on the energy extraction feasibility at Tg Tuan Headland. Goh eral. [22] further investigated the resulting 
environmental effects due to energy extraction at the Strait of Malacca, by developing the Delft3D 
hydrodynamic model 

All of the above-mentioned research studies are focused in West Malaysia only, in which the 
characterizations of tidal energy sources, as well as their feasibility analysis, are sufficient for the onsite tails 
of prototype TECs. However, not much attention is given to study the tidal steam energy resources in East 
Malaysia. Rigit er al. (11] higher tidal flow occurs atthe Sarawak in East Malaysia. Rigit er al. (11] is the 
numerical approach based on the tidal stream stations in Sarawak, in which Trio Island is identified as the 
site where the maximum tidal current speed is noted. However the mapping by {11} without a hydrodynamic 
model is not sufficient for resource characterization, Therefore for current study, Delf3D numerical 
application is utilized for setting up the hydrodynamic model representing Sarawak state, in East Malaysia, to 
characterize the tidal energy resources along the coastline of Sarawak. In the first stage, a 2D hydrodynamic 
‘model is set up to understand the depth-averaged velocity profile. Then the power density is calculated to 
find out the suitable locaton of higher kinetic energy flux. 
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4.1. Model setup 
Sarawak is the largest state of Malaysia, that is located on Borneo island in East Malaysia. It has a 

large coastline covering approximately 720 km (measured from google maps). The one side from the north is 

connected to the South China Sea, while the other sides share tertial borders with Indonesia and Burnie. 

In order to investigate the tidal hydrodynamic characteristics, Delft3D numerical tools are utilized to 
build the regional model for tidal stream energy. In the first step, Arakawa C-grid is generated in sigma 
coordinates at the model area, as represented by Figure 5. The total grid length from west to east covers 574.5 
km, and following the methodology by [9], the gird resolution close to the land boundary was set as 0.9 x 11 
km, which was gradually inereasing, and the grid resolution of 0.9x3 km is set at the open boundary. The 
flow at the area of interest is not disturbed by the instabilities at the open boundary as the distance is kept 
sufficiently far from the land boundary. After setting up the grid, the bathymetry for the proposed study area 
is derived from the general bathymetric chart for oceans (GEBCO) database for the oceans [28], 

‘The bathymetry data were interpolated to the computational grid by Delft3D-QUICKIN, which is 
shown in Figure 6. The bathymetry at the western coast of Sarawak varies around 100 m. Whereas, the deep 
sea starts from the norther-east corer of the model domain, The bathymetry at this part of the South China 
sea is up to 2000 m, this comes under the vicinity of Brunei and Sabah borders of the Borneo island. South 
China sea close to Brunei and Sabah is deeper than 2000 m, which can be a potential source for the ocean 
thermal energy conversion (OTEC). 
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Figure 6. The bathymetry of the study area and locations of observation points 
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‘The time step for the simulation was set as 1 min, as followed by the other researchers [9]. The 
‘mode! has three open boundaries named West, North, and East, as shown in Figure 6. Initial conditions at the 
‘open boundaries were set as null free water levels. Whereas the boundary conditions at the open boundaries 
are derived from the TPXO database for tidal constituents [29], 

‘The amplitude and phase of the eight major constituents M2, $2, N2. K2, K1, O1, PI, and QI are 
presented in Table 1. It can be observed that the magnitude of tidal constituents decrease from west to east 
along the Sarawak coastline. Such variations are mainly experienced due to regional bathymetry. A total of 7 
observation points are selected to observe the tidal speed at these locations along the coastline, these are 
shown with names in Figure 6. The model is run for 15 days from Ist May 2021 to 15th May 2021, to cover 
the normal neap and spring tidal hours. The final simulation findings were analyzed using QUICKPLOT 
module of Delft3D, for power density analysis and model validation, 


‘Table 1. Boundary conditions at the open boundaries of the model, a total of 8 tidal constituents are derived 


from the TPXO database 
Const Warr West, Nosh A Nook ESTA Ea 
Amplitude (an) Phase Amplitude (nn) Amplitude (my Phase Amplitude (m) Phase 
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5. RESULTS AND DISCUSSION 
5.1. Validation 

‘The infield data for the observation location in this model, shown in Figure 6, is not available. 
‘Therefore the model is validated with the water level result from the TPXO 9.2 TMD model. The TPXO 
database is used for generating the boundary condition as well as for validation of the simulation results. The 
database is established by Oregon State University, by applying Laplace Tidal formulae to the Ocean data 
acquired from the satellites, namely TOPEX/Poseidon and Jason [29]. However, TPXO is a global model, 
‘whereas in this study, we set up a regional model. The data for the two observation points named Pulau Triso 
and Bintulu can not be derived from the model, as these locations are close to the coastline which appears as 
land in the TPXO model, 

While the simulation results from the observation point Tanjung Po, Kuala Rajang, Kuala Paloh, 
‘Tanjung Sirk, and Kuala Igan are compared with the water level from the TPXO. Validity of model has been 
ensured from value of RMSE and the correlation coefficient between database results and the simulation 
results. The results are presented in Table 2, in which a good agreement between simulation and database can 
bee seen, None of the observation locations has a correlation coefficient below 0.9, while the Tanjung Sirik 
showed a very good agreement having the highest correlation coefficient of 0.96. The comparison between 
the simulation and database result of Tanjung Sirik is also shown in Figure 7. The simulation result is shown 
by the blue-colored continuous line and the TPXO result is indicated by the red-colored dotted line. 


52. 


idal flow 
‘The hydrodynamic model was impelled via density changes, local bathymetry, and tidal forces from 
the open boundaries. The current hydrodynamic model in this study is 2D, and hence the result of the 
velocity is depth-averaged. The depth-averaged velocity at Pulau Triso, Bintulu, Kuala Rajang, and Kuala 
Paloh is presented in Figure 8. While Figure 9 shows the depth-averaged velocity at Tanjung Po, Kuala Igan, 
and Tanjung Slik. From the time series analysis of both the figures, it can be concluded thatthe tidal current 
flow along the East Sarawak is very low. That's why itis not feasible to generate the power. 
Figure 8(a) shows the timeseries analysis of depth-averaged velocity (m/s) at the location Triso 
Island. The maximum depth-averaged velocity at this location is 0.25 mvs. The location Bintulu is the lowest 
site along the Eastem Sarawak coastline, and the maximum depth-averaged tidal velocity at this is up to 0.25 
rs during neap tide hours, as shown in Figure 8(b). Whereas among all observation sites, Kuala Rjanag, 
received the highest tidal current velocity of 0.55 m/s, as shown in Figure 8(c). The maximum depth- 
averaged velocity at Kuala Paloh is 0.5 m/s, as shown in Figure 8(d). Figure 9(a) shows the time series data 
of depth-averaged velocity at Tanjung Po, where maximum velocity was noted as 0.53 mis. The results of 
depth-averaged velocity at Kuala Igan and Tanjung Sirik are shown in Figure 9(b) and Figure 9(c), 
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respectively. The increase of velocity at these locations is because of the unique geographic coastline at 
‘western Sarawak, The coastline converges roughly to a single place called Triso island, as a result, the water 
is accelerated towards this site. However. the water depth at these locations is not high enough to install the 
turbines, the deepwater location is Kuala Igan. The above-mentioned locations are in the vicinity of higher 
velocity, however, there is variation in the spatial distribution of the depth-averaged velocity 


TPXO Output 


Time (Month of May) 


Figure 7. Comparison between simulation and ‘TPXO database results of water level at Tanjung Sirk. The 
red-colored dotted line is water level data extracted from TPXO TMD model 9.2 


‘The spatial analysis of tidal stream velocity is presented in Figure 10. The results are during the neap 
flood tide during the peak flow speed, as it will help to identify the locations of higher tidal flow speed. The 
‘wo main regions of higher depth-averaged velocity are identified in the vicinity of the observation points 
Tanjung Sirik and Tanjung Po. From the temporal analysis, the peak Velocities at these locations were noted 
to be up to 0.55 m/s. However, the spatial analysis shows the depth-averaged velocity near Tanjung Po is up 
to I m/s and near Tanjung Sirik is 0.8 mvs. The available tidal velocities of these locations are not suitable for 
power generation on a larger scale, but the stand-alone TECs can be installed for small-scale power 
generation 


{Pont Fae 


Figure 8. Depth averaged tidal strel velocity at the location (a) Pulau Triso, (b) Bintulu, (c) Kuala Rajang, 
and (d) Kuala Paloh 
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Figure 9. Depth averaged tidal stream velocity atthe location (a) Tanjung Po, (b) Kuala Igan, and (c) 


‘Tanjung Sivik 


‘Tabel 2. The validation results of the observation points 
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Figure 10, Spatial variation of depth-averaged tidal stream velocity in the model 
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5.3. Power flux 
‘The kinetic energy flux or the tidal power flux density gives an instantaneous idea for site selection 
atthe study area, It is calculated by the instantaneous mean Tidal power in (6). Tidal power given in (6) 


©) 


‘The factor P/A indicates the available kinetic energy flux at the study area. p is the density of the 
seawater, which was taken 1024 kg/m3, and wis the unbounded tidal stream speed. The available spatial 
kinetic energy flux (kW/m2) is shown in Figure 11. It needs to be noted that the available kinetic energy flux 
does not mean the extractable kinetic energy. The extractable kinetic energy flux depends upon the 
coefficient of performance of the tidal current turbine and the environmental factors. Nevertheless, the 
available kinetic energy flux gives the first appraisal to identify the suitable site locations. It can be observed 
that the kinetic energy flux is only shown for area within vicinity of observation locations Tanjung Sirik and 
Tanjung Po, whereas at the other locations tidal current speed is too small which is not suitable for power 
generation. The available kinetic energy flux near Tanjung Sirik ranges around 0.3 kW/m2. However, for the 
area near Tanjung Po, it ranges from 0.4 to 0.5 kWim2. Some areas with energy flux above 0.6 kWim2 can 
also be noted. Therefore, this location can be suggested for the detailed tidal energy assessment study. 
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Figure 11. The available kinetic energy flux (kW/m2) atthe study area 


6. CONCLUSION 

Oceans represent a huge potential for tidal energy harnessing. Among other renewable energy 
sources, tidal power generation has highest predictibility. The theoretical worldwide tidal current potential is, 
7800 TWhiyear. However, « lot of its potential is yet to be explored. Hydrodynamic modeling is the key 
feature in exploring the macto tidal stes. In this study, a hydrodynamic model is developed for the case study 
of Sarawak Malaysia, to explore the macro tidal sites. For this, Delft3D numerical application is utilized to 
solve the Navier Stokes equations. The model is validated against the database results, and macro tidal sites 
are identified in the study area. The two locations of Tanjung Sirik and Tanjung Po are noted to be suitable 
sites along the Sarawak coastline. The depth-averaged tidal stream velocity was noted to be up to | m/s and 
the highest available kinetic energy flux was 0.6 kW/m2. 
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